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We study a mixture of strongly interacting bosons and spinless fermions with on-site repulsion 
in a three-dimensional optical lattice. For this purpose we develop and apply a generalized DMFT 
scheme, which is exact in infinite dimensions and reliably describes the full range from weak to strong 
coupling. We restrict ourselves to half filling. For weak Bose-Fermi repulsion a supersolid forms, 
in which bosonic superfluidity coexists with charge- density wave order. For stronger interspecies 
repulsion the bosons become localized while the charge density wave order persists. The system is 
unstable against phase separation for weak repulsion among the bosons. 
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Ultracold atomic gases confined in optical lattices pro- 
vide a new laboratory for investigating quantum many- 
body problems with high precision and tunability [3, Q . 
In this way new light can be shed on notoriously diffi- 
cult problems in condensed matter physics One of 
the intriguing aspects of cold atoms is that the atomic 
quantum statistics can be controlled. In particular, cold 
atomic gases offer thepossibility to realize mixtures of 
fermions and bosons [j, H, @, 0, 0] , recently also in an 
optical lattice 0, U(3| ■ It is thus possible to create phys- 
ical systems without analog in conventional solid state 
physics. 

An exotic quantum phase that has intrigued re- 
searchers for a while is the supersolid with superfluid 
order, i.e. broken U(l) symmetry, and coexisting par- 
ticle density wave order. It is still an open question 
whether a supersolid has been realized in recent experi- 
ments on 4 He (llj . While in single-component quantum 
gases supersolids can only be stabilized by inclu ding near- 
est neighbor repulsion between the particles [l2j, they 
can be conveniently realized in Bose-Fermi mixtures with 
on-site repulsion as we show in this Letter. Earlier the- 
oretical studies already suggested that Bose-Fermi mix- 
tures can be unstable against charge density wave (CDW) 
and supersolid order or phase separation (PS). However, 
so far all theoretical approaches either dealt with one- 
dimensional systems 13 , 14 , 15 , Trj 17l. [l8l. Hij , or relied 
on weak-coupling approximations 2C , 21 A22 1 . 

Here we introduce and apply a generalized dynami- 
cal mean-field theory (GDMFT) that treats this prob- 
lem in a fully non-perturbative way. In this method the 
fermions are described by dynamical mean field theory 
(DMFT) [2^ , into which the bosons are incorporated by 
means of the static Gutzwillcr decoupling approximation 
[24| of the hopping. This approach therefore reproduces 
the strong coupling behavior of both the fermions and 
the bosons. In particular it is able to describe the for- 
mation of a bosonic Mott insulator state for strong re- 
pulsion between the bosons at integer filling. Here we 
restrict ourselves to half filling for both the bosons and 
fermions ((rib) = ( n f) = ?)■ For weak interspecies repul- 
sion we predict the formation of a supersolid phase, in 
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FIG. 1: (Color online) Schematic structure of different T — 
phases of a spinless Bose-Fermi mixture in an optical lattice. 
The red (dark) particles correspond to the fermions, while 
cyan (light) particles denote the bosons. In both the super- 
solid and the AMI phase the bosons and fermions have an 
alternating density pattern as depicted in (a). In the super- 
solid (b) the density oscillations are small and the bosons are 
superfluid. In the AMI + CDW phase (c) the density oscil- 
lations are large and the bosons are localized. The schematic 
structure of phase separation is depicted in (d). 



which the bosons form a superfluid with spatially modu- 
lated density and the fermions form a CDW. For stronger 
repulsion between the bosons and the fermions there is 
a first order phase transition to an alternating Mott in- 
sulator (AMI) plus CDW, in which the bosons become 
localized at every second lattice site while the fermionic 
charge density wave persists. For weak bosonic repulsion 
we find an instability towards phase separation (PS). We 
depict all these different phases schematically in Fig. [TJ 
A mixture of fermions and bosons in an optical lattice 
can be described by the single-band Fermi-Bosc Hubbard 
model: 
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where cj (h\) is the fermionic (bosonic) creation oper- 
ator at site i, while n{ = c\ a Ci a (n\ = b\hi) denotes 
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the number operator and the chemical potential for 
fermions (bosons). Ub and Ufb are the on-site boson- 
boson and fermion-boson interactions respectively. 
denotes summation over nearest neighbors, and t/(&) is 
the tunneling amplitude for fermions (bosons). 

Following the very succesful DMFT [23J and Gutzwillcr 
24j schemes, which are exact in infinite dimensions, we 
consider the infinite-dimensional limit (d — > oo) of this 
model first, which is expected to be a good approxima- 
tion to three spatial dimensions. In order to retain a 
finite kinetic energy the hopping parameters are rescaled 
as tf/y/d and tb/d. We then follow the standard "cav- 
ity derivation", i.e., we consider a single impurity site 
and integrate out all other sites [23| ■ In the limit of infi- 
nite dimensions the only terms that survive in the effec- 
tive action for the impurity site are the local terms, plus 
a bosonic source field and a fermionic (Weiss-) mean- 
field. The bosonic part corresponds to the Gutzwillcr 
approximation, whereas the fermionic part corresponds 
to DMFT. Therefore, the GDMFT employed in our cal- 
culation consists of the DMFT algorithm for the fermions 
[23I ] , combined with Gutzwiller mean-field theory for the 
bosons. Subsequently the action for the impurity site 
is mapped onto a generalized single impurity Anderson 
model (GSIAM). As usual, the impurity site is coupled 
to a noninteracting fermionic bath, which provides a self- 
consistent dynamical (Weiss-) mean field [23j |. In addi- 
tion, the GSIAM now also contains a bosonic degree of 
freedom, which is self-consistently coupled to the super- 
fluid order parameter, according to Gutzwiller mean-field 
theory [2i| . In summary, the GSIAM is described by the 
following Hamiltonian, which allows for a two-sublatticc 
structure: 
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Here a is the sublattice index (a = —cr), z is the lat- 
tice coordination number, (p a = (b a ) is the supcrfluid 
order parameter, and Vko- are the fermionic hybridiza- 
tion matrix elements. The hybridization function is de- 
fined as A a (uj) = it J2k I Vk CT | 2 (5(ccJ — £k)- For conve- 
nience we perform our calculations on the Bcthe lat- 
tice, which has a semielliptic non-interacting density of 
states p(e) = 2\/D 2 - e 2 /irD 2 . Here D = 2*/zt f is the 
non-interacting fermionic half-bandwidth. In the follow- 
ing we take D as the unit of energy. The fermionic 
DMFT self-consistency relation on the Bethe lattice has 
the form A a (w) = jAg-(uj) 0, where A a (uj) is the lo- 
cal fermionic interacting (impurity) spectral function. To 
calculate A a (uj) and ip a = (b a ) from the GSIAM @ wc 



use the nonperturbative numerical renormalization group 
(NRG) technique [25|, [26[. The resulting hybridization 
function as obtained via the spectral function A a (u>) and 
the bosonic order parameter tp a determine the new co- 
efficients of the GSIAM. This procedure is iterated until 
convergence is reached. The GDMFT approach incorpo- 
rates the local correlations between bosons and fermions 
in a fully non-perturbative fashion and thus reliably de- 
scribes the full range from weak to strong coupling. In 
our calculations we use a cut-off for the number of bosons 
on the impurity site, which can be kept low due to the re- 
pulsive interactions, which suppress multiple occupancy 
of the bosons. All of the results presented here are ob- 
tained at T = 0. 

The self-consistent GDMFT procedure as described 
above can yield multiple stable solutions. To find the 
ground state of the system, we need to compare the en- 
ergy of these solutions, which is given by 
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where Hk/N = | 2 CT (Hy) — ztbtfi-itpi and the indices 
±1 correspond to the two different sublattices. To calcu- 
late the fermionic part of the kinetic energy we use the 
same approach as for an antiferromagnetic state, which 



also has a two-sublattice structure [23l. \27\ 
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where p(e) is the fermionic non-interacting density of 



states and A(e,u>) = — — 
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is the spec- 



tral function, with Q a = lo + /if — T. a (uj) (cr = 
±1). We calculate the self-energy following [28| : 
= UfbF(ui)/G(u>) where G(u>) is the fermionic 
Green's function and F(uj) = {jb^bf^)^. 

We now first apply our GDMFT procedure to the limit 
Ub = 00, i.e. hard-core bosons. We fix the number of 
bosons and fermions at half filling ((rib) = ( n f) = 5)' 
which makes the system particle-hole symmetric. With- 
out loss of generality, calculations are performed for re- 
pulsive Fermi-Bose interactions: Ufb > 0. The case of 
attractive interactions will be inferred later on with the 
help of a (staggered) particle-hole transformation. Since 
we take the non-interacting fermionic half-bandwidth D 
as the unit of energy, the bosonic hopping amplitude tb 
and the interaction Ufb are the remaining adjustable pa- 
rameters. 

Our results are shown in the Ufb — tb phase diagram 
in Fig. [2] For weak repulsion between fermions and 
bosons we obtain a supersolid phase with a small CDW 
amplitude(see Fig. [3|). For strong interactions between 
fermions and bosons we obtain a bosonic alternating Mott 
insulator phase (AMI) together with a charge density 
wave (CDW) of the fermions. In this phase the fermionic 
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FIG. 2: (Color online) Phase diagram of the Fermi-Bose Hub- 
bard model with spinless fermions and hard-core bosons at 
half filling. We identify the supersolid phase (below red solid 
line), the alternating Mott insulator (AMI) phase with charge 
density wave (CDW) (above red line), and the coexistence 
region (between the red and blue (dashed) line). Energies 
are expressed in units of the non-interacting fermionic half- 
bandwidth D. 




FIG. 3: (Color online) Amplitude of the CDW for fermions 
(blue circles, solid line) and hard-core bosons (red crosses, 
dashed line) as a function of the fermion-boson interaction U fb 
for the case when ztt = 0.4D. In the inset we plot the bosonic 
supcrfluid order parameter as a function of the fermion-boson 
interaction [7/6. 



CDW amplitude |A7V/| is almost maximal, while the 
bosons are completely localized and have a CDW am- 
plitude equal to |AiVf,| = 0.5. Taking into account vir- 
tual bosonic particle-hole excitations beyond Gutzwillcr 
would however lead to a slightly smaller bosonic CDW 
amplitude. This transition is very similar to the one 
for bosons in a superlattice: upon increasing the poten- 
tial difference between the sublattices there is a Mott- 
insulator transition at half filling (29l |. For intermediate 
coupling both solutions are stable within GDMFT. To de- 
termine which of them corresponds to the ground state, 
we have compared their energies as given by Eq. ([3]). 
We find that the supersolid phase always has the lower 
energy, i.e. the ground state is the supersolid. The coex- 
istence of GDMFT solutions is a strong indication for a 
first order phase transition (at T = 0). As shown in Fig. 
[2J the critical value Wj h for the phase transition from the 
supersolid into the AMI phase increases with the bosonic 
tunneling amplitude. 

The fermionic spectrum is always gapped. Spectral 
densities are shown in Fig. [4j The gap is small for the 
supersolid phase, but at the transition point there is a 
jump in the gap and in the AMI phase it becomes of the 
order of the non-interacting half-bandwidth D (see inset 
of Fig. [3}. This implies that the latter phase will be 
more stable against finite temperature effects. 

So far we have considered repulsive interactions 
between bosons and fermions. To see what hap- 
pens for attractive interactions Ufb < we apply a 
staggered particle- hole transformation to the fermions, 
Ci — > (— l)*cj, which leads to a minus-sign in front of the 
Bose-Fermi interaction term. This implies that for at- 
tractive interactions we obtain the same quantum phases, 
but the CDW-oscillations arc now in-phasc, instead of 
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FIG. 4: (Color online) Fermionic spectral function in a mix- 
ture with hard-core bosons and zti — OAD for different values 
of Ufb- In the inset we plot the size of the gap in units of D 
as a function of Ufb- The gap is defined by the frequencies 
for which the spectral function has half its maximal value. 



out-of-phase as for repulsive interactions. 

We now proceed by considering finite interactions be- 
tween the bosons, i.e. relaxing the hard-core condition, 
but still assume the fermions and the boson to be half- 
filled. We consider the case that the bosons are slightly 
slower than the fermions: zt b = OAD. Our findings are 
summarized in the Ufb — U b phase diagram in Fig. [5] For 
strong bosonic repulsion Ub the results are similar to the 
ones found for hard-core bosons: we find a supersolid for 
weak U fb and the alternating Mott insulator phase for 
stronger Ufb, separated by a first order transition. The 
critical interspecies repulsion at the transition between 
supersolid and the AMI phase increases when the value 
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FIG. 5: (Color online) Phase diagram of the Fermi-Bose Hub- 
bard model with spinless fermions. Both fermions and bosons 
are half-filled and ztb = OAD. Stable phases are the super- 
solid (left of the red solid line and above the green dash-dotted 
line) and the alternating Mott insulator (AMI) phase with 
charge density wave (CDW) (right of the red solid line and 
above the green dash-dotted line). In the area between the 
red solid line and blue dashed line both solutions are stable. 
Below the green dash-dotted line phase separation (PS) takes 
place. 



of the bosonic repulsion Ub is reduced. This is because 
the supcrsolid state acquires a lower energy when Ub is 
decreased, whereas the energy of the AMI phase remains 
the same. For weak interactions Ub among the bosons, 
the half-filled state is unstable towards PS. In this pa- 
rameter regime we do not find a converged GDMFT so- 
lution where the bosons and the fermions arc half-filled. 
To establish the occurrence of phase separation we also 
performed calculations away from half filling. We found 
a pronounced jump in the density as a function of the 
chemical potential and coexisting solutions close to the 
position of the jump. Moreover, we observed that for 
strong interspecies repulsion the phase separation is al- 
ways complete. This allowed us to compare the energies 
of the PS- and AMI states, which yields the green line 
as depicted in Fig. [5J We have checked that comparison 
of energies yields the same boundary for phase separa- 
tion as deduced from the disappearance of a converged 
homogeneous GDMFT solution. 

Also in this case we can infer the effect of attrac- 
tive Bose-Fermi interactions by performing a staggered 
particle-hole transformation for the fermions. Phase sep- 
aration turns then into phase separation of bosons and 
fermionic holes, which is equivalent to clustering of the 
bosonic and fermionic particles. So for weak repulsion Ub 
among the bosons a system with attractive interspecies 
interaction U fb will maximize its density in part of the 
system, leaving the rest unoccupied. 

In conclusion, we have studied a mixture of half-filled 
spinless fermions and bosons in a three-dimensional opti- 
cal lattice at zero temperature. We established the pres- 
ence of a supcrsolid at weak Bose-Fermi repulsion. For 



strong interspecies interaction a first order phase transi- 
tion occurs towards a state where the bosons are localized 
and form an alternating Mott insulator. An instability 
towards phase separation was observed for weak interac- 
tion among the bosons. 
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